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Recently discovered exact integrability of zero-dimensional replica field theories [E. Kanzieper, 
Phys. Rev. Lett. 89, 250201 (2002)] is examined in the context of Ginibre Unitary Ensemble of 
non-Hermitean random matrices (GinUE). In particular, various nonperturbative fermionic replica 
partition functions for this random matrix model are shown to belong to a positive, semi-infinite 
Toda Lattice Hierarchy which, upon its Painleve reduction, yields exact expressions for the mean 
level density and the density-density correlation function in both bulk of the complex spectrum and 
near its edges. Comparison is made with an approximate treatment of non-Hermitean disordered 
Hamiltonians based on the 'replica symmetry breaking' ansatz. A difference between our replica 
approach and a framework exploiting the replica limit of an infinite (supersymmetric) Toda Lattice 
equation is also discussed. 



I. INTRODUCTION 



How replicas arise. — In physics of disorder, all observables depend in highly nonlinear fashion on a stochastic 
Hamiltonian hereby making calculation of their ensemble averages very difficult. To determine the latter in an 
interactionless system, one has to know spectral statistical properties of a single particle Hamiltonian Ti. contained 
in the mean product of resolvents, G{e) — tr{e — 'H)~^. Each of the resolvents can exactly be represented as a 
ratio of two integrals running over an auxiliary vector field ip which may consist of either commuting (bosonic) or 
anticommuting (fermionic) entries. In the random matrix theory limit, when a system Hamiltonian is modelled by an 
N X N random matrix Ti of certain symmetries, the resolvent G{e) equals 

G(e) = i J Vi^ViPiji'iPie-'^^^'''''''''^ l^j PV^PV e"'^"'^' '^''^'^ (1) 

where Sn — 'tpg {edge' — Ti-ie') ^'t, is an iV-component vector ip'^ = {ipi, • • • , ^Pn), "0 is its proper conjugate, and 
Ime 7^ 0. Summation over repeated Latin indices is assumed. 

Although exact, this representation is a bit too inconvenient for a nonperturbative averaging due to the awkward 
random denominator. To get rid of it, Edwards and Anderson^ proposed a replica method based on the identity 

_ 2 

In Z = lim . (2) 

Upon assigning to Z a meaning of a quantum partition function 

Z{e) =1^1 VipViP e-'^-^^'-'^^-'^l (3) 



the average resolvent {G{s)) can be determined through the limiting procedure 



(G(e)) = hm i I- (Z"(£)) (4) 



involving the average of the partition function 

|n| 



•' a=l 



* Published in: Frontiers in Field Ttneory, edited by O. Kovras, Ch. 3, pp. 25 - 54 (Nova Science Publishers Inc., New York, 2004). ISBN: 

1-59454-127-2. 
t Electronic address: eugene.kanzieper@weizmann.ac.il 



2 



describing n identical noninteracting copies, or replicas, of the initial disordered system The nature of replicated 
fields '(/''■"•' in ^ is determined by the sign of n in they are bosonic for n < and are fermionic otherwise. Clearly, 
the procedure assumes mutual commutativity of the replica limit n 0, an ensemble averaging (•••), and the 
differentiation operation d/de. 

Contrary to the representations and Q contain no random denominator hereby making a nonperturbative in 
disorder calculation of the resolvent G{e) viable. Depending on the origin of auxiliary fields ip^"'^ in ©, such a disorder 
averaging supplemented by identification of physically relevant low lying modes of the theory, would eventually result 
in effective replica field theory^ (also called a nonlinear replica a model) defined on either a noncompact'^''* or a 
compact^ manifold. 

Why the replica limit is problematic. — Seemingly innocent at first glance, the above field theoretic construction 
appears to be counterintuitive and rising fundamental mathematical questionsSi. Indeed, due to a particular integration 
measure which makes no sense for n other than integers, the average of (jSJ cannot directly be used to implement the 
replica limit determined by the behaviour of Z„(e) = (Z"(e)) in a close vicinity of n = 0. To circumvent this 
difficulty (which reflects^ a true, continuous geometry^ of replica field theories), one may at first evaluate the average 
replica partition function Zn{s) for all n € Z+ (or Z~), and then analytically continue the result into a vicinity of 
n = in order to make taking the replica limit well defined and safe. This route, however, is full of pitfalls. 

In the context of mesoscopic physics, the subtleties involved in carrying out the replica limit were discussed at 
length in Refs^*2iifliiiiiSii^. All these studies have debated whether or not the nonperturbative^'^ sector of replica field 
theories is reliable. This issue is of conceptual importance yet is not pure academic because the replica field theories 
are among a very few means available to address problems involving both disorder and interactions, about which the 
famous Efetov's supersymmetry approach (SUSY) - a prime tool in studying noninteracting disordered systems for 
the last two decades - has nothing to say beyond a weak- interaction limi^ . 

Approximate treatment of replicas. — In the early studjiS by Verbaarschot and Zirnbauer, a nonperturbative 
sector of nonlinear replica a models was thoroughly examined in the context of the Random Matrix TheorjiiS. (RMT) . 
Having mapped the problem of eigenvalue correlations in the Gaussian Unitary Ensemble (CUE) of large random 
matrices onto both bosonic and fermionic replica field theories, these authors had found that the two formulations 
of nonlinear replica cr model supplied different results for the density-density correlation function, both apparently 
differing from the correct one firmly established by other methodsiSiii. These findings led the authors to conclude 
that the replica method is 'mathematically ill founded'. The failure of the replica method to correctly account 
for all nonperturbative contributions to a physical observable was attributed to a nonuniqueness of the analytic 
continuation of replica partition functions in the replica parameter n away from (either negative or positive) integers. 
This standpoint, recently reiterated by Zirnbauer^^, has formed a prevailing opinion in the literature that the replica 
method may at best be considered as a perturbative tool not being able to reproduce truly nonperturbative results 
accessible by alternative SUSY technique. 

The paper that challenged the opinion about inner deficiency of replica field theories and triggered their further 
reassessment was that of Kamenev and Mezard^. Based on ideas of replica symmetry breaking originally devised in 
the theory of spin glasses^^, these authors came up with a procedure that eventually produced nonperturbative results 
for the CUE density-density correlation function out of fermionic replicas, albeit in an asymptotic region describing 
evolution of a quantum system at times not exceeding the Heisenberg time. (Subsequently, this approach was applied 
to a number of problems such as the energy level statistics in disordered metallic grains beyond^° the RMT limit, 
spatial correlations in Calogero-Sutherland modeXs^^^ , a microscopic spectral density of the Euclidean QCD Dirac 
operator—, and energy level fluctuations in Ginibre ensembles of non-Hermitean random matrices^^) . 

Briefly summarised (for a detailed exposition the reader is referred to original publicationsSiiS as well as to a 
critical analysia^^ by Zirnbauer) , the framework^ rests on an approximate saddle point evaluation of replica partition 
functions represented in terms of |n|-fold integrals containing a large parameter. In doing so, nontrivial saddle 
point configurations with so-called broken replica symmetry have to be taken into account in order to reproduce 
nonperturbative results. While leading to asymptotically correct expressions for spectral fluctuations in the Gaussian 
ensembles possessing unitary, orthogonal and symplectic Dyson's symmetries, the procedureSi cannot be considered as 
mathematically satisfactory because it involves a nonexistingi^J^ analytic continuation of replica partition functions 
to a vicinity of n = 0, the domain which is crucially important for implementing the replica limit. 

Towards exact integrability of replica cr models. — Analysis of Refs&SiiSiiiiiS (see also RefiM) hints that 
approximate evaluation of replica partition functions is the key point^— to blame for inconsistencies encountered in 
the procedure of analytic continuation away from n integers. In such a situation, leaning towards exact calculational 
schemes in replica field theories is a natural move. 

A step in this direction was taken in the recent paper—, where partition functions for a number of fermionic replica 
a models were shown to belong to a positive, semi-infinite Toda Lattice Hierarchy extensively studied in the theory of 
nonlinear integrable systems24. In conjunction with the r-function theory^^'^®-^^'^^ of the six Painleve transcendents 
(which are yet another fundamental object in the theory of integrable hierarchies), this observation led to exact 
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evaluation of replica partition functions for a number of random matrix ensembles in terms of Painleve transcendents. 
Resulting nonperturbative Painleve representations (which implicitly encode all hierarchical inter-relations between 
partition functions with various replica indices) were used to build a continuation of Z^s away from n e Z+. While 
not addressing the important issue of uniqueness of such an analytic continuation, the route of Ref.^'^ has yielded - 
for the first time - exact nonperturbative results for random matrix spectral statistics out of fcrmionic replicas. 

More recently, Splittorff and Verbaarschofe22i2ii have suggested that such nonperturbative results could directly 
be obtained from the replica limit of an infinite Toda Lattice equation without Painleve reduction whatsoever. In 
fact, the approachlMS, rests on the observation that, if properly normalised, the fermionic and the bosonic replica 
partition functions of a zero-dimensional interactionless system form a single, infinite (that is, supersymmetric) Toda 
Lattice Hierarchy belonging to its either positive (fermionic) or negative (bosonic) branch. While greatly simplifying 
calculations of spectral correlation functions through a remarkable bosonic-fermionic factorisation^^, the framework 
developed by Splittorff and Verbaarschot is, to a large extent, supersymmetric in nature as it explicitly injects^^ a 
missing bosonic (or fermionic) information to otherwise fermionic (or bosonic) like treatment. 

In the present paper, a detailed account is offered of a nonperturbative approach^ to zero dimensional fcrmionic 
replica field theories which is based on exact Painleve representation of replica partition functions with the emphasis 
strongly placed on technical details. Specifically, we focus on the Ginibre ensemble'^^ of complex random matrices with 
no further symmetries. This particular random matrix model is of special interest in the light of recent findings that 
associate statistical models of normal random matrices with integrable structures of conformal maps and interface 
dynamics at both classical^'^ and quantum scales'^^. (The reader is referred to Rcf.''^ for an introductory exposition 
of these recent developments and to Refi^ for a review of other physical applications and extended bibliography) . 

The paper is organised as follows. In Sec. II, we collect the basic definitions and present the major results 
regarding the Ginibre Unitary Ensemble of random matrices. In Sec. Ill, a fermionic replica field theory approach 
to non-Hermitean complex random matrices is outlined and integrability of the field theory is established. The 
integrability which manifests itself in emergence of a positive, semi-infinite Toda Lattice equation for replica partition 
functions and also results in exact representations of replica partition functions in terms of Painleve transcendents, 
eventually culminates in reproducing exact fluctuation formulas for Ginibre complex random matrices. Finally, in 
Sec. IV, we make a comparison of our exact approach with the approximate treatment^^ of non-Hermitean disordered 
Hamiltonians based on the 'replica symmetry breaking' ansatz; we also comment on differences between our approach 
to replicas and a framework^^''^" exploiting a replica limit of the infinite Toda Lattice equation. Most lengthy technical 
details are collected in the two appendices, A and B. 



II. GINIBRE UNITARY ENSEMBLE OF NON-HERMITEAN RANDOM MATRICES: DEFINITIONS 

AND BASIC RESULTS 

Preliminaries. — Statistical ensemble of generic N x N complex random matrices H £ C^x^ whose entries are 
independently distributed in accordance with the Gaussian law2i 

Pjv(H) = TT-^' exp (-tr HH^) (6) 

has first been introduced in the pioneering work by Ginibrc"^^. Throughout the paper, such random Hamiltonians will 
be denoted as H £ GinUE^r. 

The joint probability distribution function Pn{zi, • • • , zn) of N eigenvalues of the matrix Ti. is of particular interest. 
Having Schur-decomposed the Ti. a.s Ti, = W {Z + TZ) U where t7 is a unitary matrix which is unique up to the phase 
of each column, 7?, is a strictly upper-triangular complex matrix, and 2^ is a diagonal matrix Z = diag (zi, • • • , zn) 
consisting of N complex eigenvalues {zi} = {xi + iyi\ of the TL, Ginibre managed to derive the joint probability 
distribution function oi {zt} in the form 

N N 

Pn{zi,--- ,zn)^C{N) II \z,,-zeJ^Y[w^{ze,ze) (7) 

where the weight function w'^{z,z) equals w^(z,z) — exp(— zz). Given the integration measure (PZi = dxgdyi, the 
inverse normalisation constant in iQ is determined to be C~^(iV) — 'k^T{N + 1). 
Of primary interest is the p-point correlation function 

ATI r ^ 

i?p(zi,...,Zp;7V) = — — — / W (fZ.P^iz^,--- ,zn) (8) 
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describing a probability density to find p complex eigenvalues around each of the points zi , • ■ ■ , Zp while positions of 
the remaining levels stay unobserved. Quite often, one is also interested in the thermodynamic limit of the correlation 
function 

Pp(zi,--- ,Zp) = lim i-Rp ( T^>--- :T^;^ ) (9) 

that magnifies spectrum resolution on the appropriate energy scale while letting the matrix size N tend to infinity. 

The multi-fold integral in © can explicitly be evaluated by adopting the Gaudin-Mehta-'^^ method of orthogonal 
polynomials originally introduced in the context of Hermitean random matrix theory. It is a straightforward exercise 
to demonstrate that i?p(zi, • • • , Zp; N) admits the determinant representation 

Rp{zi, ■■■ ,Zp;N)^ det [K]y{zk, Zi)]^.^^^^^,,, _p (10) 

involving the scalar kernel 

Kn{z, z') = w{z, z) w{z', z') Pt{z)PfXz') 

expressed in terms of polynomials Peiz) orthonormal in the complex plane z — x + iy 

^ d?Zw^{z,z)Pk{z)Pt{z)^5ut (11) 

with respect to the measure i(j^(z, z) cP Z . 

For instance, it follows from (I1U|) that the density of states and the two-point correlation function equal 

Ri(z-N) = Kn{z,z) (12) 

and 

R2{zi,Z2]N) = Kn{zi,Zi)Kn{z2,Z2) - \Kiq{zi,Z2)\^ , (13) 

respectively. 

For the Gaussian measure, the orthonormal polynomials Pi{z) are just monomials 

P,{z) - ^ (14) 

V^r(£ + i) 

leading to the scalar kernel 

^A^(^,^') = ie-/2,-.'.72^ iItT)- (15) 
Density of states. — Put into the integral form, the kernel H15|l yields the finite-A^ density of states 

-22 /-oo 

i?i(z;7V) = — — / dXe-\\ + z-zr-\ (16) 
In terms of the upper incomplete gamma function 

/>oo 

T{a,x) = / dtf'-^e-' (17) 

J X 

the level density equivalently reads 

A careful analysis of the integral H16|) shows that, in the large-iV limit, N complex eigenvalues are (almost) uniformly 
distributed within a circle of the radius \/N centered at z = 0, 



i?i(2;Af > 1) ~7r-i6'(V^- |z|), (19) 
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d{x) being a Heaviside step function. 

In the vicinity Zc — {VN + u) e^'^ of the edge \z\ = \/N of the two-dimensional eigenvalue support described by 
(|19|l . the density of states sharply crosses over from Ri{z;N ^ 1) = tt^^ at \z\ < \/N to Ri{z;N 3> 1) = at 

\z\ > \fN . The crossover is described by the local density of eigenvalues i? J*^''^"* (u) = Ri{zc] N) which, in the large- 
limit, turns out to be independent of the matrix size N , 

1 

4"'^'(«) = ^^ / dte-'''\ (20) 
Expressed in terms of the complementary error function 

2 2 

erfc(a;) = ^ / dte"*' (21) 
V"" Jx 

the local density of states reads 

- ^crfc (u\/2) . (22) 
At |u| ^ 1, the tail asymptotics can be read out of H22|l and are given by 

«<••■■•'(.) 4 + 5^)^ (23, 

Two-point correlation function. — At finite N , this simplest fluctuation characteristics is given by (|13|l and H15|l . 

In the large- limit, when both zi and Z2 in (|13|l are situated inside the circle \z\ < ^/N, the kernel (|15() reduces to 
simple exponentials so that the two-point correlation function becomes A'^-independent, 

i?2(^i,^2) =i?2(^i,^2;A^» 1) = ^ (l-e-l^i-^^l') . (24) 



Make notice that the two-point correlation function R2{zi, Z2; N) differs from the density- density correlation function 
R{zi, Z2; N) defined by 

R{z,,Z2;N) = {tr6'{z^-n) tr d\z2 - H))^^^.^^^^^ . (25) 
The latter contains an additional (5-function contribution^ 

R{zi, Z2) = R{zu Z2\ N)^- 5'^{zx - Z2) + i?2(zi, Z2) (26) 

TT 

coming from the self-correlation of eigenlevels that chanced to meet in the complex plane. 



III. FERMIONIC REPLICA FIELD THEORY 



In this section we are going to re-derive the above nonperturbative results by mapping the Ginibre ensemble © 
of complex non-Hermitean random matrices onto a zero-dimensional (OD) fermionic replica field theory. Since the 
nonperturbative fluctuation formulae collected in Sec. ^1 ^^'^ weW known for almost four decades (we remind that 
Ginibre's work'^^ dates back to 1965), one may wonder why we should bother ourselves with such a minor issue. 
The answer prompted by the discussion in Sec. ^is twofold. First, more than twenty years after their invention^, 
replica field theories largely remain unexplored territory from the viewpoint of their controllable treatment away from 
a perturbative sector. Second, learning intrinsic integrable structure of fermionic replica field theories in the simplest 
OD limit - apart from encountering indisputable mathematical beauty of exact theory - creates a basis for future work 
beyond the RMT: Given an intimate connectionSS, between integrability and the underlying physical symmetries, one 
may hope that some crucial characteristics of OD replica partition functions are in fact not so peculiar to these simple 
models but remain true in a more general setting. 
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A. Replica partition functions 

Density of states. — To determine the average density of complex eigenvalues of the matrix Hamiltonian H e 
GinUEjv, we use a proper modification^^ of Q and jSjl- A new recipe^ is 

i?i(z;7V) = lim — -^Z„(z,z;Ar) (27) 

n^Q 7771 OZOZ 

where the replica partition function Zn{z, z; N) is determined by the matrix integral 

Z„gK(z, z; N) = (dct "(z - n) det "(z - H^))^^^,^^^^ ■ (28) 

Representing each of the determinants in (|28|1 as a field integral over an A'^-component fermionic field and perform- 
ing the averaging over Ti. e GinUETv, one maps Zn{z,z;N) onto a fermionic replica sigma model, Zn{z, z; N) 
Zn{z,z-N), of the form^ 



(z,z;7V) = (det^( ) ) ■ (29) 



z 



QeGinU„ 

Importantly, while Zn{z, z; N) in H28|l is defined for arbitrary n G M, the replica parameter n in the representation H29(l 
for the mapped replica partition function Zn{z, z; N) is restricted - by derivation - to positive integers only, n G 
To emphasise this difference between two types of partition functions, we will write either 'tilded' Z„ {n G Z+) or 
'untilded' Z„ (n G R). 

The matrix integral H29|) can be reduced to an n-fold integral^^ by making use of a singular value decomposition of a 
complex matrix Q G C"^". Expressing it as Q = UKV where U G U(n)/U(l)", V G U(n) and A = <iia.g{\\''^ , • • • , A,V^) 
with \i > 0, and calculating a Jacobian of the transformation Q — > (A, C/, V), one derivesSS 

/>oo ^ ^ 

Z^{z,z;N)= l[dXee-^^{Xe + zzf J] |A,, -A,,^. (30) 

•^0 1=1 £l<£2 = l 

We reiterate that this representation makes sense for n G Z+ only. This is precisely the reason why the replica limit 
(|27|l with Zn replaced by Z„ cannot be implemented directly. 

Density- density correlation function. — Similarly to the level density (|27|l . the density-density correlation hmc- 
tion H25fl can be retrieved from the replica limit 



R{zi, Z2\N) = lirn^ a,, az: n-, a;: Zn{zi,Zi]Z2, Z2] N) (31) 



ri^o tt'^tt? dzidzi dz2dz2 
where, for arbitrary n G K, the replica partition function Z„(zi, zi; Z2, Z2; N) is defined by 

Znmizu ~zi\ Z2, Z2\ N) = (det "(zi - H) det "(zi - H^) det "(z2 - H) det "(z2 - ^^))«eGinUE„ ■ (^2) 

Upon a fermionic mapping^ which restricts the replica parameter to n G Z+ , this generating function can be rewritten 

as 

Znez+{zi,zi;z2,Z2;N) ^ (dct'^ [ ^ ]) (33) 



QeGinUE2„ 



with the diagonal matrix Z = diag(zi 1„, Z2 1„). Note a tilde in (|33|l . 

For zi, Z2 finite and of order unity, Znez+ can be reduced'^^ to a matrix integral over U G U(2n) which eventually 
boils down to the n-fold integral^ 



n^^A,e(-^/2)A, -Q \Xe,-XiJ'. (34) 

- 1 i> 1 n ^ I) 1 



^1 1=1 ei<£2=i 



Here, 

Zl + Z2 



LO = Zi — Z2- (35) 



As is the case (|3(J|) . this representation makes sense for n G 2+ so that the replica limit (|31|l with Z„ replaced by Zn 
cannot be implemented directly. 
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B. Replica partition functions as members of a positive Toda Lattice Hierarchy 

By derivation, the ?i-fold integral representations (|30() and H34|l of the rephca partition functions Zn{z,z;N) and 
Zn{zi, zi; Z2, Z2', N) stay vahd for n G Z+ only. Therefore, any attempt to retrieve spectral fluctuation properties of 
the matrix Hamiltonian Ti. out of the replica limits H27|l and H31|l with Z„ replaced by will inevitably face the 
problem of analytic continuation^*^ of Z^s away from n positive integers. For this procedure to be controlled, an exact 
result for Z„ is desired to start with. 

For approximate treatment^'^- of the above replica partition functions the reader is referred to Ref*^. In this 
subsection we wish to explore another route which rests on exact and, therefore, truly nonperturbative evaluation of 
replica partition functions. A connection between nonlinear replica a models and the theory of integrable hierarchies 
is at the heart of our formalismS^. A proof that the nonperturbative fermionic replica partition functions form a 
positive, semi-infinite Toda Lattice Hierarchy is the first important outcome of our approach. 

Bulk density of states. — To show how the Toda Lattice Hierarchy emerges in the context of (|30|) . we represent 
the Vandermondc determinant there as 

n 

n (A,, ~A,J =det(A^,-')M=i,-,«, (36) 
ei<i2=i 

simultaneously shift all Af's therein by zz, and perform the n-fold integral (|3U|I by means of the Andreief-de Bruijn 
integration formulc^i*^ 



/ f[dn{Xe) det[Ak{Xe)ht=i^... ,ndet[Bk{Xt)]k,e=i,- ,n = n\ det ( f dfi{X)Ak{X) Be{X)) 



(37) 



which holds for any benign integration measure dfi{X) given convergence of the integrals involved. Up to (irrelevant 
for our purposes) factorial prefactor, the interim result is 

Zn{z,z;N) = det ( f dXe~^{X + zz)^+''+^') . (38) 

V"'0 J k,l=0,--- ,n-l 

While exhibiting some beauty (in particular, e~'^^ Zi(z, z] N) coincides, up to a normalisation prefactor, with the 
density of states i?i(z; + 1) in GinUEjv+i, see (fTH|l l the representation is not very informative or helpful. 
What is more helpful is another though totally equivalent form of (|38(l . 

Z^{z, z; N) = e"^^"(zz)"("+^) f„(^z; A), (39) 

which involves the Hankel determinant 

f„ (zz; A) = det [a*;+^n (zz; A)! (40) 

with 

fo(zz;A) = l, n(zz-; A) = j'^ dX A^e^ = ^^^^-j^+r^ ■ (41) 



Here, T{a,x) is the upper incomplete gamma function (|17|l . The initial condition for fo(zz; A) reflects the fact that 
Zo(z, z; A) = Zo(z, z; A) = 1, see (EHJ. 

The Hankel determinant H40() is a remarkable object. Whatever the function fi (zz; A) is, by virtue of the Darboux 
Theorem^, the entire sequence {r^jg^+l satisfies the equation 



^1 ^2 



r. 



fn-lTn+U n&'L+ (42) 



where the prime ' stands for 9(^2)- Equations H39() and (|42() taken together with the known initial conditions f o = 1 
and fi given by l|41|l establish^^ a hierarchy between nonperturbative fermionic replica partition functions Zn with 
different n G Ij^ . The exact result H39I) , (|41() , and \A2\ is an alternative to an approximate positive- integer-n treatment 
of the very same replica partition function presented in Ref.--. 

Equation {4°/^ , known as positive, semi-infinite Toda Lattice equation in the theory of integrable hierarchies^, is 
the first indication of exact solvability hidden in replica field theories. Importantly, emergence of the Toda Lattice 
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Hierarchy is eventually due to the (3 = 2 Dyson's symmetry of the fermionic replica field theory encoded into the 
squared Vandermonde determinant in l|3U|) . 

Density of states at the edge. — The above symmetry argument ensures that the replica partition function for the 
edge density of states in the GinUEjv will obey the same Toda Lattice equation albeit with different initial conditions. 
Close to the edge \z\ = \/N (see discussion next to (|19|l '). one is interested in the large-iV replica partition function 

Z„(z, z; N) taken at z = {\fN + u) e'"^ in the regime \u\ <C \fN . The latter, denoted as Zn'^^^^\u), is independent of 
the matrix size N and equals 

Zt''^'>iu)= J|dife-*'/2-2ut, -Q \ti,-ti,f. (43) 

Irrelevant numeric prefactors were omitted. 

Much in line with previous calculations, the n-fold integral H43|) can be transformed into the Hankel determinant 
form 



Z^^'^^'^Xu) = det \dt+'z['^'''\u)\ _ ^ (44) 

with 



fc.£=0,--- ,n-l 



= 1, Zf'''''\u) = e^"' j dte-'"/^ = ^1 e2"'erfc (uV2) . (45) 

As soon as Zq^^^^'' (u) = 1 (normalisation), the Darboux Theorem^ can be applied to conclude that the entire sequence 

} of replica partition functions at the edge of the two-dimensional eigenvalue support belongs to the positive, 
semi-infinite Toda Lattice Hierarchy 

^(tails) ^(tails) _ (^(tails) ^ ^(taUs) ^(taUs) ^ ^ ^+ ^^g^ 

the prime ' stands for du- Also, similarly to the observation made below (|38|l . we notice that, up to a prefactor, 
g-2M2^(taiis) j.^^ coincides with the edge density of states r'^'^^^^\u) as given by (|20|l . We will comment on this later 
on. 

Bulk eigenvalue correlations. — By the same token, the replica partition function (|34() designed to calculate the 
density-density correlation function via the replica limit H31|l belongs to a Toda Lattice Hierarchy, too. Separating z 
and w-dependent pieces in one derives 

^n(^i,^i;;^2,^2;A^) = e~^"^^~''^ Z,,{iuQ) (47) 
where Zn{uJoj) determined by the n-fold integral in (I34II can be cast into the Hankel determinant formi^iS^ 

Z„(.^)=detf5|=:4^^^^i^) . (48) 



LULj/2 



k,e=o,--- ,71-1 



The positive, semi-infinite Toda Lattice equation for {^fegz+} readily follows by virtue of the Darboux Theorem. 



C. Replica partition functions and Painleve transcendents 



While important for revealing integrability of the field theory, the Toda Lattice equation for fermionic replica parti- 
tion functions .^„gz+ is not much helpful in performing the replica limit, if taken alone. Indeed, a positive Toda Lattice 
equation gives no close expression for Z„ as a function of n e Z+ that would facilitate an analytic continuation to 
n G M in general, and to the region < n <C 1 in particular. As was recently shown by Splittorff and Verbaarscho1iS2i^, 
this difficulty can be circumvented if one succeeds in gaining a complementary hierarchical information about bosonic 
replica partition function. (The latter satisfies a negative? Toda Lattice equation). Viable yet surprisingly efficient^ii 
for random matrix models describing interactionless stochastic systems, this route is certainly unavailable for exact 
replica description of physical systems in presence of interaction^ which requires to use either bosonic or fermionic 
field integrals in order to properly accommodate quantum statistics of interacting species. 

Considering exact replica treatment of disordered interacting systems as a legitimate goal, it would be conceptually 
important to not rely on such a complementary information. Fortunately, for OD interactionless systems at hand. 
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it is indeed possible. Miraculously, the same Toda lattice equation governs the behaviour of so-called r-functions 
arising in the Hamiltonian formulation2SiSS42LS& of the six Painleve transcendents (PI - PVI), which are yet another 
fundamental object in the theory of nonlinear integrable systems. Luckily, the Painleve equations being second order 
nonlinear differential equations contain the hierarchy (or replica) index n as a parameter. As will be demonstrated 
below, this feature of Painleve equations makes them serve as a proper starting point for constructing a consistent 
analytic continuation of nonpcrturbative replica partition functions away from n integers. This Painleve reduction 
further confirms exact solvability of replica a models and assists^^ performing the replica limit. 

Bulk density of states. — Certainly being an option, the aforementioned Toda to Painleve reduction'*''^ is not the 
only way to arrive at the sought Painleve representation of the replica partition function Zn{z, z; N). An alternative 
approach would rest on the observation that the n-fold integral (|3U|I is essentially a Fredholm determinant'^^ associated 
with a gap formation probability 

within the interval (0, zz) in the spectrum of an auxiliary n x n Laguerre unitary ensemble. Celebrated result*^ due 
to Tracy and Widom states that 



El°--)(0;a) = exp(^'^rfi^) 



(50) 



where av{t) = cr„(t; a) is the fifth Painleve transcendent satisfying the Jimbo-Miwa-Okamoto form of the Painleve V 
equationSSii£ 

{ta'^)^ - (adV)^ - ((TV - ta'y) \av - ta'y + 4a'y (^a'y + n+ = (51) 
supplemented by the boundary conditio n^^i^° 

2 

'^vit)\t^+oo--rit + an-^^ + 0{t-^). (52) 
We, thus, derive an exact Painleve V representation of the replica partition function Zn{z, z; N) in the form 



zz 



Z^{z, z; N) = e"^^" exp ( / dt ) (53) 



t 

where crv(t) = an{t;a — N). 

Note that (|51|l ~ (|53|l contain the replica index n e Z+ as a parameter. Taken together with the fact that the above 
Painleve representation encodes all hierarchical inter-relations between the replica partition functions with various 
replica indices, it is very tempting to conjecture that H53|) . as it stands, holds beyond n e Z+ as well. Indeed, in 
Appendices A and B, it is proven that (|51|l - H53() stay valid for generic real valued n > — 1 so that the replica limit 
(|27|l with Zn substituted by Zn can safely be implemented. 

To proceed, we expand a solution to (|51|l around n = 0. Owing to the normalisation Zq = I, the expansion starts 
with a term linear in n, 

oo 

ay{t) = a„(i; a) - ^ nPfp{t; a). (54) 
p=i 

Only the first term of the above series is of our interest since the replica limit relates the bulk density of states 
Ri{z; N) to the function fi{t;N) as 

Riiz;N)^7T-'[l + f[izz;N)]. (55) 
Here fi{t;N) satisfies the differential equation 

{tf^ f - (/i - tf[r - 2Nf[{f, - tf[) - {Nf[ f = (56) 
subject to the conservation constraint 

poo 

d^zRi{z;N)^ dt[l + f[{t;N)]^ N. (57) 
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Identifying a complete square in (|56|l . we reduce the latter to the Kummer differential equation 

h + {N~ t)f[ ± tf'l = 0. (58) 

The constraint H57I) makes us look for those solutions fi{t;N) whose first derivative /{ is bounded at t = +oo and 
possibly has an integrable singularity at t = +0. This class of functions sought welcomes the sign (— ) in (|58|l leading 
to a general solution 

hit;N) = C\iN -t) + C2t''+' iFiiN,N + 2;-t). (59) 
The constraint H57I) uniquely fixes unknown constants Ci and C2 be Ci = and C2 = — l/T{N + 2). This yields^ 

/i(^;^)-- p(^^^) iFi(jV;jV + 2;-t) (60) 

where iFi{a;b;t) is the confluent hypergeometric function of Kummer. Consequently, the replica limit H27|l of the 
fermionic partition function (|53|l results in the density of states (I55|l of the form 

which is equivalent to H16|) and (|18|l . The small-n expansion of the fermionic replica partition function 



In Zn{z, z; N) = n{zz) 



[z-z) 



N 



{N + l)T{N + 2) 



2F2{N,N + 1;N + 2,N + 2: -zz) 



(62) 



is behind the result H61|) . It should be stressed that the finite- result (|61|l cannot be produced by approximate 
treatmentiSS, of replicas which heavily relies on availability of a large parameter (N ^ 1) in the integral representation 
(EOJl. 

Density of states at the edge. — With the replica partition function Z„ given by (|43|l . the density of states 
at the edge \z\ = ^/N of the two-dimensional eigenvalue support is determined by the replica limit 



i?r^)(u)=lim^|^Zr'^)(u). 



(63) 



The partition function zi''*''^'' can again be viewed as a Fredholm determinant associated with a gap formation 
probability 



/u n n 



(64) 



<^2 = l 



within the interval (w, 00) in the spectrum of an auxiliary n x n Gaussian Unitary Ensemble. In terms of the fourth 
Painleve transcendent^^ ctiv, it reada^SiS 



EKoo)(o) =E(f^°°)(0) 



exp 



dt CTiv (t) 



where (Tiv(i) = cr„(i; a — 0) satisfies the Painleve IV equation in the Jimbo-Miwa-Okamoto form 

((tJ'v)' - Hta[y ~ aiv)' + Mv{<^'w ' '^a){a[y + 2n) - 
subject to the boundary condition^ 

n{a + n) 



<^iv(t)lt 



-2nt 



(65) 



(66) 



(67) 



In both (|66|l and H67() the parameter a has to be set to zero. The above equations result in the Painleve IV represen- 
tation of the fermionic replica partition function 



Zi*''"''(M) = 



e^"" exp 



nV2 



dtaiY{t) 



(68) 
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which holds for n G Z"*". 

To perform the repHca Hmit we follow the technology that led us to the small-n expansion (|62|1 . To this end 
we have to assume^ that the Painleve IV representation l|()8|l stays valid in a vicinity of n = 0. Writing down^ 

oo 

criv(t) = o-„(t; a) = ^ n^gpit; a), (69) 
p=i 

one derives from here, H()3|) and (|68|l that 

zL r 1 _ 1 

(70) 



TT 



l + ^g[{-uV2; 0) 



Equation for 51 = gi{t; a ~ 0) follows from (|66|l and H()9|) . 

± 2(tg'i - gi) = (71) 

while the boundary conditions are^ 

9ih°o) + 2 = (convergence) 
g[{+oo) + 3/2 — (conservation) ^ ' 

The sign (— ) in (|71|l leads to a solution with unbounded first derivative g'i{t) at both infinities and is therefore 
incompatible with H72() . Equation H71(l with the sign (+) yields a general solution 

c,i(t) = Cit + C2 ( ierft+ ^e-*' ) . (73) 



Boundary conditions (|72|) fix the constants^ be Ci = — 7/4 and C2 = 1/4. By virtue of H70|l this results in the density 
of states 

4*^"^)(^) = -Lerfcfw%/2V (74) 
This is identically equivalent to H23|l . The small-n expansion of the fcrmionic replica partition function 



lnZ(*^"^'(^^) = - 

is behind the result (|74(l . Being exact, the formula H74|l describes the tails of level density both inside (u < 0) and 
outside (u > 0) of the circle^^ \z\ = VN- The approximate treatment^! of replicas has failed to reproduce the density 
of states outside the circle, \z\ > \/N . 

Density- density correlation function. To determine this spectral characteristics, we put the replica limit l|31|) 
into the form 

^(-i'-2)= lim Qaf-a^) i^^dl-di) t„(z,z-;^,^) (76) 

involving the variables z and lo as defined by (|35|l . The notation T„(z,z; 0-1,(1)) stands for the replica partition 
functional. 

'fn{z,z\uj,u)) = Z„ + (77) 

At z ± ijj/2 of order unity, the n-fold integral representation (|34(l makes it possible to express T„ in terms of a gap 
formation probability 

E(r^'-)(0;a)= / Y[dXee-^'X1 Y[ |A,, -A,,p (78) 

i = l il<i2 = l 



/27r 



: ue 



erf 



+ 0(n2) 



(75) 
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within the interval (ww, oo) in the spectrum of an auxiliary n x n Laguerre unitary ensemble (compare to (|49|l ). The 
result due to Tracy and Widom'*^ states that 

where the fifth Painleve transcendent (Ty(t) = an{t;a = 0) satisfies the equation (|51fl with a — and meets the 
boundary condition^ 

-v(i)|,_-^e-*(l + Ori)). (80) 

The above equations yield, for ri G Z+, the exact representation 

. - -^ exp [ni2zz + u;Lu/2)] ( r,.ayit)\ 
T„(z,z;w,u;) = , exp - / — ^ . (81) 

To implement the replica limit, we have to analytically continue H81() into a vicinity of n = 0. Although at the 
moment we do not have a proof that (|81|l as it stands also holds for n away from positive integers, armed with the 
previous experience we are going to conjecture that this is indeed the case so that 

Ri^u - ^ + 4 li- 4 id.d^f ^^^^^ exp (- r dt ^) . (82) 
Equation H82|) suggests that the two functions, hi{t) and /i2(i), of the small-n expansion 



ay{t)=J2n''hp{t) (83) 
p=i 

contribute the density-density correlation function in the replica limit. As we tend to avoid explicit reference to the 
boundary conditions for Painleve transcendents, the easiest way to determine hi is to notice that, at \z\ <C Vn, the 
bulk density of states equals 

i?i(z;iV) = lim — -^Z„(z,z;0,0;JV). (84) 
n^o nn ozoz 

This is so because, at n ^ 0, the partition functions H28|l and 132|) taken at zi = z and Z2 = become indistinguishable 
if considered as functions of the energy variable z. Given (|77|l and (|81|) . one derives 

Z.{z. -z. 0, 0; N) . exp (- £ dt ^) . (85) 

Only linear in n term of the expansion (|83|l contributes the replica limit (|84|) yielding 

^i±Z^. (86) 

TT 

According to the replica result (|61|l . this must be equal to I/tt whence we conclude that the function h\(t) is a 
constant. Further, the equation 

(th![f - {hi - th[f ^ (87) 

following from 151|) and H83|l sets 

hi(t) = 0. (88) 
Therefore, the first nontrivial term in the expansion (|83|l actually starts with •n?h2{t) where h2(t) satisfies the equation 

th'^ ± (th'2 ~ h2) ^ 0. (89) 
The sign (— ) is the one that meets existence arguments applied to (|79|) . As a result, we come down to 

h2{t) = Cit + C2 E2{t) (90) 
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where E2 (t) is the exponential integral 



E„{z) = J dt—, Rez >0. (91) 

Again, by existence arguments, Ci must be set to zero to ensure convergence of the integral in the exponent of l(7i^ . 
To fix the constant C2, we make use of the observation that T„ has to be finite at w = 0. This brings C2 = 1 so 
that^ 

h2{t)^E2{t). (92) 

Collecting (|81|l , H83|l , l|88|l and H92|l , we end up with the following nonperturbative small-n expansion^^ of the logarithm 
of the replica partition function: 

In t„(z,z;w,il)) = n (2zz + ^) - ji^[\n{ujLj) + Ei{ujlu) - E2iujLo)] + 0{n^). (93) 
The replica limit H7()|) applied to (|93|l culminates in the exact result for the density-density correlation function 

i?(zi, Z2) = -S^izi ~ Z2) + 4 f 1 - . (94) 

Notice a presence of the (5-functional contribution in (|94|) describing the self-correlation of complex eigenlevels. The 
latter is inaccessible by the approximate treatmentSS, of replicas. 

D. Comment on a puzzle 

In the paper—, Nishigaki and Kamenev have noticed that there exists a 'striking resemblance' between the finite- 
N densities of states and the replica partition functions for non-Hermitean random matrices belonging to all three 
universality classes /3 = 1, 2 and 4. In the context of the present study, the observation of the authors of Ref— can 
be translated into the identity (see two remarks below and (^BJ, respectively) 

i?i(z;7V) (X e-^^Zi(z,z;iV- 1) (95) 

that links the density of states Ri{z; N) in GinUEAr to the replica partition function Zi{z,z;N — 1) for the same 
ensemble albeit of a smaller dimension GinUE7v_i. 

The identity (|95|l is not a miracle and can well be understood as a consequence of the relation (|12|l . Indeed, viewing 
the scalar kernel Kn{zi, Z2) in H12|) as a matrix integral^2iSS 

KM{zuZ2)^e-^--^-/^e-^^-^^/^ (det (zi - H) det (z2 - Wt))„^^.^yj,^_^ 

one readily identifies 

Kn{z,z) (xe-^^ Zi{z,z;N -I) (96) 

whence (|95|) foUows^^. Interestingly, so explained identity (|95|l taken together with H27|l leads, in the context of 
fermionic replica field theory, to a much less trivial statement 

—— Z^{z, z; N) = ™ e-""Zi(z, z;N~l) + ©(n^). (97) 
ozoz 

Similar small-n expansions should exist for two other {13 — \ and 4) universality classes in non-Hermitean RMT. 



IV. DISCUSSION 



In the present paper we have offered a detailed account of a nonperturbative approach^'^ to zero-dimensional 
fermionic replica field theories which is based on exact representation of replica partition functions in terms of Painleve 
transcendents. Focussing on Ginibre ensemble'^^ of complex non-Hermitean random matrices, we have revealed an 
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intrinsic integrability of associated replica field theories. It materialises in two ways: First, at n G Z"*", the replica 
partition functions were proven to belong to a positive, semi-infinite Toda Lattice Hierarchy. Second, the very same 
replica partition functions were shown to be expressible in terms of solutions to Painleve equations which (i) contain 
the replica index as a single parameter and which (ii) implicitly encode all hierarchical inter-relations between the 
fermionic replica partition functions with various replica indices. 

The above two observations [(i) and (ii)] led us to conjecture that Painleve representations of fermionic replica 
partition functions stay valid beyond n e Z+ and, in particular, in a vicinity of n = 0. Indeed, for a particular 
case of the replica partition function Zn{z, z; N) designed to determine the bulk density of complex eigenvalues, this 
conjecture was rigorously proven by appealing to the Okamoto r-function theory of the Painleve V. (Similar, in 
spirit, proof was previously given^^ in the context of the one-point Green function in the finite-iV Gaussian Unitary 
Ensemble where a Painleve IV equation arises). Once justified, it is no surprise that taking the replica limit of the 
Painleve-represented fermionic partition function Zn{z, z; N) has culminated in reproducing exact nonperturbative 
results for the bulk density of states. In other cases, which include the tails of level density and the density-density 
correlation function in the spectrum bulk, although implemented without a formal justification, the replica limit 
of replica partition functions expressed in terms of Painleve transcendents has also brought exact nonperturbative 
results. This fact as well as other encouraging applications^ of the present method make us look further into the 
rational reasons behind its success. Possibly, recent developments^ in the field of extended Toda Hierarchy may give 
us the right lead. 

Finally, a remark is in order aimed to pinpoint the difference between the approach^ detailed in this paper and a 
complementary approach recently elaborated in RefsiSS*^. Although both approaches exploit the very same replica 
representations of quantum correlation functions as a starting point, the two frameworks are conceptually different. 
The present approach ''■^'^ based on exact Painleve evaluation of fermionic replica partition functions followed by their 
continuation into a vicinity of n = makes no reference whatsoever to bosonic partition functions. On the contrary, the 
approach22i^ exploiting the replica limit of the Toda Lattice equation for replica partition functions rests explicitly, 
and unavoidably, on the observation that suitably normalised fermionic and bosonic replica partition functions are 
the members of a single Toda Lattice Hierarchy albeit belonging to its positive (fermionic) and negative (bosonic) 
branches, respectively. Implemented on the level of such an infinite - supersymmetric in essence — Toda Lattice 
equation, the replica limit reveals a remarkable factorisation of quantum correlation functions for an interactionless 
matrix Hamiltonian into a product of both fermionic and bosonic partition functions. It is this factorisation which - in 
order to be materialised on the operational level - explicitly infuseg^i^S^ a missing bosonic (or fermionic) information to 
what early appeared to be a pure fermionic (or bosonic) formulation of the field theory. While facilitating calculation 
of quantum correlation functions in the interactionless case, this feature makes the approacbSSiffi be potentially 
inapplicable for an exact nonperturbative replica treatment^ of OD Hamiltonians with interaction8SMS4S& whose 
presence requires using of either bosonic or fermionic field integrals in order to accommodate a proper quantum 
statistics of interacting species. 
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APPENDIX A: REPLICA PARTITION FUNCTION Zn{z,z) FOR THE DENSITY OF STATES BEYOND 

n e z+ 

To prove that the Painleve representation H53|l . H5U|) . (|51|l and 1)52(1 stays valid for generic real- valued n > —1, we 
start with the definition of the replica partition function l|28|l . Upon Schur decomposition of the matrix Ti G GinUEjv 
specified below the partition function H28|l can be put into the form 



Acknowledgments 




(Al) 



15 



Having introduced a new set of integration variables = — z, we use (|36|) and a proper extension of the Andreief-de 
Bruijn formula (|37|) into the complex plane to derive 

^„(z, z; N) =det( f Si C""^ T"^^ e-(5+^)(€+2)'\ . (A2) 

\J J k,i=0,--- ,N-l 

Decomposition ^ = pe'^ helps us perform integration in (|A2(I eventually resulting, for n> —1, in 



Z„(z,z;iV) = 



'^^^ det 



dpe-P%^^+^+''+'h-i{2\z\p) 



(A3) 



Here, Ik{z) is the modified Bessel function. No care was taken of numeric prefactors turning to unity in the replica 
limit n —f Q. 

The {k^t) matrix entry in IjASp can be recognised as the confluent hypergeometric function of Kummer 



Applying the transformation formula 

iFi (a; b; z) ~ iFiib — a; b; — z) 
and performing trivial row and column operations under the sign of determinant, one concludes that 



Z„(z,z;iV) = det 



k,e=0,--- ,N-1 



where 



(A4) 



(A5) 



(A6) 



(A7) 



is the generalised Laguerre function. 

Importantly, the representation (|A6p holds for generic real- valued n > —1. Moreover, as is proven in the Appendix 
B, the determinant ljA6|l is identically equivalent to another one 



Z^^)(z,z;7V) = det 



L 



(k) 

n+l-k 



k,l=Q.,--- ,N-1 

arising in the r-function theorjsSSii^ of the fifth Painlcve transcendent so that 

Z„(z,z;7V) = z(^)(z,z;A^). 
In the notation of Forrester and Witte^^ (see their Sec. 3.5) the latter can be expressed as 

Zl:\z, z; N) = e-^^^" (zz)-^'/^ f[iV](zz) 
where f[N](t) is the Hankel determinant 

f[N]{t) = det [6^,+' (e* iFi(-n; 1; -t))], , 



,N-l 



(A8) 

(A9) 
(AlO) 
(All) 



built on the operator 5t = td/dt. 

In accordance with the Okamoto r-function theorjiSSii^ of the Painleve V, the f[N]{t) (which also depends on n) 
can be related to the fifth Painleve transcendent a^ 



CTv(i) = t 



dt 



In 



■f[N]{t) 



(A12) 



Here, fTv(0 — crn{t\ N) obeys^ the Painleve V equation H51() with a = N subject to the boundary condition H52|l . 
This fact taken together with the equations (jA10|) and ljA12p proves that the representation (|51|) - H53|) . as it stands, 
stays valid for real valued n > — 1 in general and, in a vicinity of n = 0, in particular. This justifies taking the replica 
limit in the way described between (|54|) and Ht)l|l . 
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APPENDIX B: A PROOF OF EQUIVALENCE OF THE TWO DETERMINANTS dXell AND lIXsll 



Throughout the Appendix, t denotes t — zz. 
(i) First, we use the identity 



to reduce (|A6|I to 



Zn{z,z;N) = det 



k,e=0,--- ,N-l 



(ii) Second, we use the identity 

_d_ 

dw 



r(^ + m+l) x-^ l'^^-^Hw) 

r(m + l) ^ ^ ^+"' 



to bring (|B2p to 



Z„(z,z;iV) =det 



k,e=0,--- ,N-1 



(iii) Third, we use the identity 



to rewrite (IB4II as 



(Bl) 



(B2) 



(B3) 



(B4) 



(B5) 



or, equivalently, as 



Z„{z,z;N) ^ det 



Zn{z,z;N) = det 



k,e=0,--- ,N-1 



ki-t) 



k,e=0,--- .N-l 



(B6) 



(B7) 



The transformations (i) to (iii) were aimed at expicit eUmination of deep cancellations occurring in the determinant 
(jA6|) as can be recognised by comparing the entries of matrices under the sign of determinant in (jA6|) and (|B7p . 
The determinant 1B7|I is already very close to the desired form (|A8|I . the only difference lying in the upper index 



of the Laguerre functions {L 
determinant, we will constantly use the identity 



{k~e) 

n+i-k 



vs To demonstrate that this difference does not affect the value of the 



(B8) 



when performing a set of row operations under the sign of determinant in (jB7|l as detailed below. 

As the first row in HB7() and the first row in HA8|I are already equivalent to each other, we start with the last, the 
A^-th row in (jB7|) . Adding to it the {N — l)-th row raises the upper index of the A^-th row entries by 1. In the next 
step, we add the content of the {N — 2)-th row to the {N — l)-th row; this also raises the upper index of the entries 
in the {N — l)-th row by 1. We go on with this procedure up unless we arrive at the row jj 2 to which we add the 
content of the row jj 1. This procedure (consisting of A — 1 steps) has 'repaired' the row tt 2 bringing it to the form 
of the row jj 2 appearing in (jA8|l . 

To proceed further, we come down to the last, A-th, row in HB7I) adding to it the content of the previous, (A — l)-th 
row. Again, we go on with this procedure up until we reach the row jj 4 to which we add the content of the row 3. 
This makes the row ft 3 fit the same row in ljA8p . This took {N — 2) steps. 

Now we again come down to the very last row and repeat the same procedure. In (A — 3) steps we will find that it 
brings the row tt 4 to the desired form. Proceeding in the same manner, we realise that A(A — l)/2 row operations 
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which do not change the value of the determinant (|B7|I bring the matrix under the sign of determinant in l|B7|l to 
that in HA8|I . This completes the proof of the formula (|A9() . 
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